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ABSTRACT.  Developed  in  this  report  are  expressions 
for  the  lift,  pitching  moment,  and  location  of  the  aero¬ 
dynamic  center  due  to  free-  stream  angle  of  attack  for 
an  isolated  cylindrical  wing  based  on  a  solution  of 
Weissinger's  general  boundary  condition  for  ring  wings. 
From  the  components  of  ring  vorticity  expanded  in  a 
Birnbaum  series  and  satisfying  the  velocity  boundary- 
flow  condition,  and  from  the  Kutta-  Joukowski  relation 
for  the  force  on  a  vortex  surface,  the  lift,  moment,  and 
aerodynamic  center  can  be  expressed  in  terms  of  the 
Birnbaum  coefficients.  It  was  found  that  the  lift  and 
pitching  moments  for  cylindrical  ring  wings  are  sub¬ 
stantially  different  from  the  results  for  two-dimensional 
flat  plate  (thin)  airfoils.  Calculations  show  a  significant 
shift  in  the  aerodynamiXjUiantajv  upstream  of  the  1/4- 
chord  point  for  cylindrical  Ling! wings'. 
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This  report  presents  a  method  for  deter¬ 
mining  the  loadings  on  an  isolated  ring  wing 
at  angle  of  attack.  The  results  of  this  method, 
based  on  linearized  lifting- surface  theory, 
were  applied  in  estimating  the  pitching  moments 
on  a  ring-wing  tail  configuration  used  as  a  con¬ 
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NOMENCLATURE 

Fourier  coefficient  of  an(^) 

Birnbaum  coefficient  of  gn(£) 

Birnbaum  coefficient  of  gj(^) 

Birnbaum  coefficient  for  unit  angle  of  attack 
Total  normal  coefficient  of  lift 
Total  pitching  moment  coefficient 
Fourier  coefficient  of  y(t,,  <j>> 

Ring  vortex  density  due  to  conicity 

Ring  vortex  density  resulting  from  cambered  section 

Ring  vortex  density  resulting  from  angle  of  attack 

Length  of  cylindrical  wing 

Total  pitching  moment 

Radius  of  cylindrical  wing 

Velocity  influence  function 

Velocity  influence  function  for  aj 

Free-stream  velocity 

Free-stream  radial  component  of  velocity 
Aerodynamic  center  measured  fi  om  mid-chord  point 
Coordinates  of  cylindrical  wing 
Pitch  plane 

Total  net  force  normal  to  wing  axis 
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a(£.  4>) 


Qn ( £ )  >  n  =  0,  1,2  ... 

°0*«) 

a0 


al 

y(i.  4>) 


Local  angle  of  incidence 
F ou.-iei  coefficient  of  a(C,  <j>) 

Slope  of  wing  profile  lor  cambered  section 
Half- cone  angle  (conicity)  of  wing  section 
F ree  stream  angle  of  attack 
Total  ring  vortex  density 


Ap 

i  =  ZX/L 
V 

r\  =  \(i  -  £') 


Force  per  unit  area 

Dimensionless  axial  coordinate  - 1  <  t  <  1 
Dimensionless  variable  of  integration  -1  <  C,1  <  ] 
Argument  of  velocity  influence  function 


9  Angular  coordinate  in  Birnbaum  expansion 

k  =  L/2R  Wing  length- diameter  ratio 

p  Fluid  density 


« 
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INTRODUCTION 


Two-dimensional  airfoil  theory  is  known  to  be  in  a  highly developed 
state  and  to  provide  an  accurate  description  of  airfoil  flow  characteris¬ 
tics.  In  ordinary  two-dimensional  thin  airfoil  theory,  two  classes  of 
problems  exist:'  (1)  the  direct  problem  m  which  the  velocity  and  pres¬ 
sure  fields  are  known  and  the  airfoil  geometry  is  to  be  determined;  and 
(2)  the  indirect  problem  in  which  the  geometry  is  given  and  ihe  flow 
field  characteristics  about  the  airfoil  are  to  be  determined.  A  typical 
problem  in  which,  the  indirect  method  Is  employed  is  the  calculation  of 
.-.paawise  loadings  on  finite  wings  for  which  the  geometry  of  the  wing 
is  known. 

Until  recently,  three-dimensional  ring  airfoil  •'heory  had  been  less 
well  developed.  However,  techniques  employing  lUting- surface  theory 
now  account,  lor  three-dimensional  effects  due  to  corvatime  of  the  chord 
plane.  It  has  teen  shown  that,  if  the  radius  of  curvature  of  the  chord 
plane  for  a  ring  wing  is  small  relative  to  the  chord  length,  the  aerody¬ 
namic  characteristics  for  a  ling  wing  differ  markedly  from  the  flow 
characteristics  of  a  two- dimensional  airfoil. 

This  report  is  concerned  with  loadings  on  the  isolated  cylindrical 
ring  wing  at  angle  of  attack.  A  subsequent  report  will  discuss  loadings 
on  ring  wings  of  arbitrary  profile  with  wing-body  interference  effects. 


GENERAL  FORM  OF  BOUNDARY  CONDITION 
FOR  RING  WINGS 

Using  cylindrical  coordinates  X,  R,  4>  with  the  origin  at  the  wing 
center  as  represented  in  Tig.  1,  a  thin  ring  wing  is  replaced  by  a 
cylindrical  ortex  surface  of  radius  R.  This  is  equivalent  to  a  thin 
airfoil-type  approximation  in  which  the  velocities  and  pressures  on  the 
ring  wing  are  represented  by  the  velocities  and  pressures  induced  by  a 
cylindrical  vortex  surface  at  the  radial  distance  R,  where  R  is  inter¬ 
preted  to  mean  the  average  radial  distance  of  the  wing  camber  line. 

The  net  loading  or  difference  in  pressure  per  unit  length  per  unit  span 
induced  by  the  vortex  surface  is  represented  by  the  Kutta- Joukowski 
relation 


Ap  =  pVy(£,  4>R  C  =  — • -1  <  l  <  1 

L 


(i) 
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FIG.  1.  Coordinates. 

where  y(£,,  41/  is  the  density  of  the  ring  vortices  on  the  cylindrical  sur¬ 
face.  The  density  can  be  written 

Y(C.  4>)  =  V  n^0  gn{l)  cos  ni>  (2) 

gn(£)  =  C0nct8  “  +  v?i  cvn  sin  v  0,  =  -  cos  0,  n  =  C,  1,  2  .  .  .  (3) 

The  components  of  ring  density,  gn(C)>  are  represented  as  a  3irn- 
baum  series,  Eq.  3,  and  thus  satisfy  the  Kutta  condition  at  the  trailing 
edge  y(l ,  4>)  =  0. 

Since  the  present  concern  is  with  the  loadings  on  the  ring  wing  for 
a  specified  angle  of  attack  between  the  wing  axis  and  the  free  stream 
velocity  and  a  specified  local  slope  of  the  wing  profile  relative  to  the 
wing  axis,  the  indirect  approach  is  taken,  from  which  the  ring  vortex 
density  distribution  is  solved  from  the  known  geometry.  Combining  the 
slope  of  the  profile  with  angle  of  attack  gives  an  effective  geometry  or 
local  angle  of  incidence 

a{£,  4>)  =  n?0  Qn(U  cos  n<t>  (4) 

an(£)  =-—  +  apn  cos  p  0,  £  =  -cos0  (5) 
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The  ring  vortex  density  distribution,  y(t,,  <j>),  must  satisfy  the 
boundary  condition  by  removing  the  radial  wash  or  local  angle  of  inci¬ 
dence,  a(f,,  4>).  For  a  fixed  value  of  n,  the  boundary  condition  in 
Weissinger's  notation  is  written 


n\  -  1 

an(£)  cos  n<j>  =  — 

4  J-  1 


gn{V)  cos  n 4>  dt, ' 


'  1  gn(£’l 

-  cos  n cp  d£' 

-l  &  -  v 


w  tie  re. 


X  rl 

+  —  gn(^')  Un(o)  cos  n  ‘P  d^'  n  =  0, 

2ir  J-l 


L 

x  =  — ,  r,  =  \a  -  v) 

2R 


2  ,  .  . 


(6) 


The  total  boundary  condition  is  obtained  by  summing  the  boundary 
conditions  for  corresponding  components  of  a (£,,  <j>)  and  y (£,,  <j>)  over  all 
values  of  n.  Thus 


nX.  f*  1 

n_0  cn(U  cos  n  9  =  —  gn(t')  cos  n  cj>  d£' 

4  o_l 

1  Z  f  1  gn^'^  cos  n  ^  d^' 

+  2t r  n  =  °  J-l  i  -  V 

\  v  f  1 

*  — nt0  gn^')Unh)cos  n4>  d^'  (?) 
2tt  1  J  -  1 

For  a  symmetrical  ring  wing,  the  local  slope  of  the  profile  is  in¬ 
dependent  of  4>.  Thus  the  local  slope  of  the  profile  is  the  term  of 
a(t,  <P)  on  the  left  side  of  Eq.  7  for  the  case  n  =  0.  Setting  ao(£)  =  ao**^) 
then,  for  a  wing  with  conicity  and  camber,  the  slope  of  the  profile  rela¬ 
tive  to  the  wing  axis  can  be  written 

a0**(U  =  t0  +  aC^>  <8) 

where  o.q  is  the  half-cone  angle  and  <1q*(£)  is  the  slope  of  the  wing  pro¬ 
file  for  the  cambered  section. 

From  Fig.  2  it  can  be  seen  that  for  any  wing  section,  the  free- 
stream  radial  component  of  velocity  is  dependent  on  <}>  and  determines 
the  remaining  term  of  a{£,,  <j>)  on  the  left  side  of  Eq.  7  for  the  case  n=  1. 
The  free- stream  radial  component  of  velocity  becomes 
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FIG.  2.  Free-Stream  Angle  of  Attack  as  a 
Function  of  4>., 


Vr  =  V  sin  a.}  cos  <j>  (9) 

For  small  cq,  Eq.  9  can  be  written 


-  =  cq  cos  4>,  a j  constant  (10) 

V 

where  cq  is  the  angle  of  attack  between  the  wing  axis  and  the  free 
stream  in  the  XZ  plane  and  cq  cos  $  is  the  angle  of  attack  in  the  XR 
plane. 

Combining  the  slope  of  the  profile  with  free- stream  angle  of  attack, 
Eq.  8  and  10,  and  substituting  in  Eq.  7,  the  total  boundary  condition  be¬ 
comes 

**  1  *  f1  ... 

Q0  **(£)■*■  Q!  cos  4>  =  |  — - _+_  g0  a')  U0(q)dC- 

2tt  «/  -  1  &  -  2m  J-  1 

\  ri  l  r  l  gi(C')  cos<|)  d£- 

+  -  gi(^')  cos  4>  dt,'  +  —  - 

4  2tt  J-l  £  - 

\  r  l 

+  —  j  gi(C')  Uj  (q)  cos  <i>  dt'  (11) 
2tt  j-  1 

where  gQ**(U  =  g0*(U  +  g0(^>  and  an^>  =  gn<?)  =  °>  n  =  2,  3,.  .  . 
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Referring  to  Eq.;  3,  5,  and  6,  the  coefficients  of  gi(£)  and  go**(£)  can 
easily  be  solved  for  in  terms  of  the  coefficients  of  cq(£)  and  qq**(?,)  from 
Tables  1  and  Z  for  the  case  n  =  0,  n  =  1.  With 

Y(C,  4>)  =  V[  g0**(£)  4-  gl(f,)  cos  4>]  (1Z) 

and  using  Eq.  1,  the  net  loading  per  unit  area  on  the  wing  surface  be¬ 
comes 

AP  =  pV2[g0**(U  +  gi(&)cos<t>]  (13) 

TABLE  1.  Matrix  Needed  To  Solve  for  CpQ  in 
Eq.  Z8  for  n  =  0 


0 

0 

-z 

0 

0 

0 

1.  080 

0 

-0.  080 

0.  5 

0.  17Z 

-Z.  187 

-0.  013 

0.  0Z7 

0 

-Z.  036 

-0.  001 

0.  009 

0..  000 

1.  ZOZ 

0 

-0.  ZOZ 

1.  0 

0.  484 

-Z.  5Z1 

-0.  08Z 

0..  118 

0 

-Z.  155 

-0.  OOZ 

0.  039 

0.  000 

0 

0.0Z8 

0 

-z.  oiz 

0 

0.  119 
0 

-Z.  054 


TABLE  Z.  Matrix  Needed  To  Solve  for  Cp^  la 
Eq.  Z8  for  n  =  1 


\ 

Mi“ 

1 

1 

0 

0 

0 

0 

-Z 

0 

0 

0.  0 

0 

0 

-z 

0 

0 

11 

0 

0 

-Z 

0.  549 

0.  381 

0.  068 

0.  003 

-0. 119 

-  1.  851 

-0.  015 

-0.  016 

0.  5 

-0.  009 

-0.  006 

-1.  98Z 

0.  000 

0.  001 

-0.  004 

0.;  000 

-1.  995 

0.  398 

0.  44Z 

0.  141 

0.  0Z1 

-0..  ZC8 

-  1.  645 

-0.  073 

-0. 076 

1.  0 

-0.  034 

-0.  038 

-1,  904 

-0.  OOZ 

0.  000 

-0.  OZZ 

0.  000 

-  1. 970 
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Recalling  that  a,j**(£)  is  the  slope  of  the  wing  profile  relative  to  the 
wing  axis,  then,  from  Fig.  1  the  normal  force  per  unit  area  can  be 
written 


dY 

-  =  Ap  cos  a (14) 

dX 

dY 

-  =  Py2  [gn**(?>)  +  gi(£)  cos  4>]  cos  a0**(£)  (15) 

dX 

dY  L 

-  =  pV2 — [gn**(U  +  gi(C)  cos  <j>]  cos  aQ**(^)  (16) 

dt  2 

where 

2X  -L  L 

i  =  ,  -  <  X  <  — ,  -1  <  £<  1 

L  2  2 


Equation  16  can  also  be  interpreted  as  the  local  lift  normal  to  the 
wing  axis  per  unit  length  per  unit  span,  and  for  an  element  of  wing  span 
Rdcf),  the  local  lift  in  the  pitch  plane  per  unit  length  (Fig.  1)  becomes 

dY 

-  cos  <(>  Rd4>  (17) 

dC 


The  total  lift  normal  to  the  wing  axis  per  unit  length  is  then 


rZir  dY 

Jo  d£ 


cos  4>  Rd4> 


f2ir  L 

=  p V2— [gQ**(t,)  +  gi(t)cos  b]cos  a**(.?,)cos  4>  Rd4>  (18) 

J  0  2 

and  the  total  lift  normal  to  the  wing  axis  beumes 

cos  <)>  Rd4>  d£ 


Y  = 


r  1  r  2ir  dY 

J-iJo 

tic 

pV2RL 

r1  ( 

2  J 

-J 

p  v2rl 

2  J 

-l  J 

p v ‘■RL  r  l  rZir 

=  -  [g0'V'*(C)+  gi(^)cos<j>]  cos  a0**(^)cos4>  d<j>  d£ 

2  J-IJO 


[gu**(^)cos<t>  cos  e.g**(£) 

+  gjJU  cos2  4>  cos  a0**(£)]  d(j)  d£  (19) 
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For  small  values  of  a q**  (Z,) 
pV2RL  C  1  CZtt 

Y  =  -  [go**(£)cos  4)  +  g,(^)cos24)]  d(j>dC  (20) 

2  J-l  J  0 


Since 


dM  L£,  dY 

d£  2  d£ 

the  total  moment  about  the  mid- chord  point  due  to  forces  normal  to  the 
wing  axis  becomes 

M  --  pV2R  |_j  |”  J  t,  [gQ**(?;)cos  $  +  g^Ocos2  4>]  d<}>  d£  (21) 

Inspection  of  Eq..  11  and  20  shows  that  go**(£)  corresponding  to  the 
slope  of  the  profile,  ao**(£)  =  +  aQ*(C).  yields  a  zero  net  force  nor¬ 

mal  to  the  wing  axis,  which  is  the  case  for  a  symmetric  ring  wing  at 
zero  angle  of  attack.  It  is  apparent  then  that  the  total  or  net  normal 
force  is  due  solely  to  the  density  distribution  gl(£)  due  to  angle  of 
attack  cq. 


LIFT  FORCE  ON  CYLINDRICAL  RING  WING 


From  Eq.  20,  the  net  force  normal  to  the  wing  axis  becomes 


Y  = 


pV^RL  pi  r  2tt 


r i  r21 

1-1 J  0 


Y  = 


J 


irpV^RL 


gl(£)cos  4>  d<j>  dt, 


gl&)  d?, 


(22) 


(23) 


Substituting  Eq.  3  in  Eq.  23  for  the  case  n  =  1 
irp V2RL  fit 


f  ( Cq i  ctg  —  +  y  Cvi  sin  vo\  sin  8  d9 

J0'  2  v  =  1  ' 


=  up  V2XJR2  I  Cqi- 


it 


1 +  cos  8 


sin  8 


sin  8  d8  +  ^  ,  C„i  sin  v3  sin8  d6j  (24) 

J  0  V  =  1  / 


where  X 


2R 


7 


NAVWEPS  REPORT  8401 


TT 


b  =-  fv(0)cosp  0d0,  r|  =  M£-  £'), 

TT  J  0 

£,  -  -  cos  0 

(29) 

r  tv 

fv(0)  =  cos  v9,Un(r|)d0, 

(30) 

JO 

For  a  cylindrical  wing,  qq**{£)  =  0,  at  angle  of  attack  ai  the  local 
angle  of  incidence,  Eq.  5,  becomes 

a01  V 

Ql  =  —  +  p  =  i  ap l  cos  pe,  p  =  1,  2  .  .  ,  (31) 

2  /•it 

aDi  =—  qi  cos  pOdG  (32) 

ir  J  0 

From  Eq.  32,  the  Fourier  coefficients  for  qj,  constant,  are  then 

apl  =  °>  P  =  1.  2  .. ... .  (33) 

a01  =  2al 


In  matrix  notation  Eq.  28  becomes 


r 

a0n 

C0n 

a  In 

=  Mn 

C  In 

apn 

^pn 

where  the  elements  of  are  calculated  from  Eq.  28  with  elements  of 
the  inverse  Mjj'l  shown  in  Tables  1  and  2  for  n  =  0,  n  =  1. 

Setting  n  =  1  in  the  matrix  equation  and  multiplying  both  sides  by 
the  solution  for  the  coefficients  Cpj,  p  =  0,1,2  ...  takes  the  form 


Coi  ' 

C1 1 

^  Mj"1 

’  aoi 

all 

.CRl. 

.  apl  . 
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From  Eq.,  33 


C01 

C11 

r— ( 

1 

jH 

S 

II 

'  2cq  * 

0 

-Si  _ 

O  * 

1 _ 

(34) 


Inspection  of  Eq.  27  and  setting  cq  =  1  in  Eq.  34  show  that  the  nor¬ 
mal  coefficient  of  lift  of  a  cylindrical  ring  wing  per  unit  angle  of  attack, 
cq,  can  be  written 


TT 

cL=- 

2 

p  (0 )  ,  1  p- 

C01  +  2  Cl] 

(0) 

c  (0) 
U01 

i—H 

1 

*— H 

2 

II 

1— 

2 

cn(0) 

0 

C2,<°> 

0 

,Cpl«»_ 

• 

where  Mq-^  is  determined  by  the  value  of 


(35) 


(36) 


Taking  X  =  1  and  recalling  a  rule  for  matrix  multiplication,  the 
values  in  the  first  column  of  Table  2  multiplied  by  the  factor  2  in  Eq.  36 
gives  a  solution  to  the  coefficients  Cpl(0),  p  =  0,  1,  2,  ...  .  Thus  for 
\  =  1,  CqiW  =  0.796,  Cqi(O)  =  -0.415,  and  for  a  unit  angle  of  attack  in 
radians,  =  0.92.  The  solution  to  Eq.  36  can  also  be  found  in  Table  3. 
Figure  3  shows  the  trend  of  the  slope  of  the  lift  curve  in  degrees  for 
different  values  of  \. 


For  the  case  of  a  ring  wing  with  X.  =  0,  Cq\(^)  =  2,  Cn(0)  =  0  and 
from  Eq.  35  for  unit  angle  of  attack  in  radians,  Cl  =  tt. 
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ABLE  3.  The  Birnbaum  Coefficients  Cpp°'  of  the 


C  (0 ) 
C01 


C,1<°> 

C21<°> 

r  (0) 
C31 

0 

0 

0 

-0.  238 

-0. 018 

0.  002 

-0.  415 

-0.  068 

! 

0.  000 

X.  0.1  £  X  £  0  6 
X.  0.6  <  X  <  1.3 


FIG.  3.  Total  Lift  Coefficient  per  Unit  Angle 
of  Attack,  deg,  versus 
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MOMENT  FOR  CYLINDRICAL  kING  WING 
From  Eq.  21,  the  total  moment  about  the  mid- chord  point  is  written 


M  = 


/  L\2  r  1  rZ-n 

,V2R|_  £gl(£)cos^4>  d4>  dC 

\  2  /  J-l J 0 


(37) 


where  from  Fig,  1  the  moment  is  negative  clockwise. 
Integration  on  *  gives 


M  =  trpV2\2R3  j  ^  £gl(£)d£ 


(38) 


Using  the  Birnbaum  series  expansion  for  gj(£)  in  Eq.  3  and,  re¬ 
calling  that  £  =  -cos  0,  Eq.  38  becomes 

M  =  irpV*-  \2  R3J"  -  |  Cqi  ctg  —  +  j  Cvj  sin  v©j  sin  0  cos  0  d0 

=  Trp  V2  \2  r3  f  .  c01 

Jo  [ 

=  -rrpV2\2R3|  - 

J  0 


1  +  cos0 ' 


sin0 


|  sin  0  cos  0  +  Cv^  sin  v0  sin  0  cos  0 

sin  20 


d0 


Cqi(cos0  +  cos^  0)  +  y  sin  v0 


=  TTp  V2  \2  R3 


r  tt 


sin220 


-  Cqj  cos  0  -  Cm  cos^O  -  C 


01 


"21 


d0 


IT  TT  1 

M  =  irp  V2  \2  R3  -  Ca  2  -  -  C21  - 

2  4 


(39) 

(40) 


Defining  the  dimensionless  coefficient  as 
M 

C jyj  =  - - - -  with  Cj^  negative  clockwise 

2ttRL2jPV2 


then  from  Eq.  40  the  total  moment  coefficient  due  to  normal  forces  can 
be  written 


Cm  =  - 


(41) 


12 


NAVWEPS  REPORT  8401 


then  from  Eq.  27  and  41,  the  aerodynamic  center  of  lift  for  ring  wings 
becomes 


result 


X  JL 

ac  _  16 

(  2Cq !  +  C21) 

TT 

L 

2 

(C01  +  iCll) 

Xac  _  i  ( 

2C01  +  C2 1 ) 

L  8  ( 

C01  +  icn) 

=  0, 

and  from  Table  2, 

together  with 

=  0. 

Equation  42  then  gives  the  well- 

L 

L 

xac  ~  ’ 

- <  X  <■ 

4 

2 

=  1, 

from  Table  2,  Eq.  34  and  42, 

(42) 


L 

2 


(43) 


C01  =  0.  398  (2aj) 

Cn  =  -0.  208  (2a2) 

C21  =  -0.  034  (2aj ) 
xac  1  (1.  592-0.  068)(4 
L  8  (0.796-0.208)^ 

Xac  =  -0.  325 L  (44) 


The  dependency  of  the  wing  aerodynamic  center  on  the  shape 
parameter,  \,  of  the  ring  wing  is  shown  in  Fig.  5.  Relative  to  a  two- 
dimensional  airfoil  of  unit- span  length,  a  chord-diameter  ratio  of  X  =  1 
for  a  cylindrical  wing  of  span  length  2ttR  will  cause  a  30%  shift  in  the  aero¬ 
dynamic  center  toward  the  leading  edge. 

CONCLUSIONS 

A  ring  wing  with  a  small  profile  slope  (Fig.  1)  is  replaced  by  a 
continuous  distribution  of  vortex  rings  of  constant  radius.  Within  the 
limits  of  linearized  theory,  the  velocities  and  pressures  on  the  cam¬ 
bered  ring  wing  with  conicity  are  represented  by  the  velocities  and 
pressures  on  a  cylindrical  vortex  sheet.  The  boundary  condition  on 
the  ring  wing  is  given  in  Eq.  11.  Inspection  of  Eq.  11  and  20  shows  that 
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J 

ft 

X 


FIG.  5.  Aerodynamic  Center 
Versus 


tlie  slope  of  a  symmetric  ring  wing  profile,  ao**(£).  gives  a  zero  net 
force  normal  to  the  wing  axis  and  that  the  total  normal  force  results 
from  free- stream  angle  of  attack.  Thus  as  a  first  order  approximation 
the  normal  force  on  a  ring  wing  of  arbitrary  profile  shape  is  equivalent 
to  the  normal  force  on  a  cylindrical  ring  wing. 

For  a  cylindrical  wing  at  angle  of  attack,  aq,  where  =  g**(£) 

=  0,  the  boundary  condition  is  defined  in  Eq.  28  for  n  =  1,  i.  e.,  the  re¬ 
maining  terms  containing  cq  and  g^  in  Eq.  11.  A  reduction  of  the 
boundary  condition  to  the  linear  algebraic  form  in  Eq.  28  is  given  in 
the  Appendix.  For  n  =  1,  with  coefficients  of  radial  wash,  api,  given, 
the  coefficients  of  the  Birnbaum  expansion  for  vorticity,  Cpf,  can  be 
determined  from  the  set  of  linear  algebraic  equations  shown  in  Eq.  28. 
This  gives  the  lift  and  moment  for  a  cylindrical  wing  as  defined  in 
Eq.  27  and  41  , 

Figure  6  gives  the  total  lift  coefficient,  Cl,  versus  cq  in  degrees 
for  different  values  of  \.  The  lift  and  moment  coefficients  per  unit 
angle  of  attack  in  degrees  are  plotted  for  different  values  of  the  ring 
wing  parameter  K  as  shown  in  Fig.  3  and  4.  The  location  of  the  wing 
aerodynamic  center  measured  from  the  mid- chord  point  is  shown  in 
Fig.  5  for  different  values  of  \.  For  X.  =  0,  the  aerodynamic  center  is 
located  at  the  1/4-chord  point,  which  is  the  two-dimensional  result.  An 
example  of  the  three-dimensional  effect  due  to  curvature  of  the  chord 
plane  is  demonstrated  for  the  case  \  =  1  which  indicates  about  a  30% 
shift  in  the  aerodynamic  center  upstream  from  the  i/4- chord  point. 
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“i 


FIG..  6.  Total  Lift  Coefficient  Versus 
aj,  deg,  for  Various  Values  of 

Tables  1  and  2,  obtained  from  E  28,  give  the  Birnbaum  coeffi¬ 
cients  required  to  calculate  the  lift,  moment,  and  aerodynamic  center 
for  isolated  cylindrical  wings.  Table  3  was  obtained  from  Eq  36  and 
Table  2.  Recent  investigations  show  a  significant  change  in  aerody¬ 
namic  characteristics  of  the  ring  wing  in  the  presence  of  a  central 
body.  It  has  been  found  that  because  of  the  presence  of  axial  forces 
further  changes  are  apparent  in  wing  loadings  and  location  of  the  aero¬ 
dynamic  center  of  lift  due  to  wing  camber  and  conicity.  The  results  of 
this  study  will  be  published  in  the  near  future. 
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Appendix 

SOLUTION  TO  THE  BOUNDARY  CONDITION 


The  general  boundary  condition  is  written 


n\  rl  i  /»l  gnU')  a.  rL 

an(U=—  gn(t')di'+—  -  d£’+—  gn(^')Unh)d^  (A-l) 

4  J-l  2-n  J-l  t,  -  Ztt  J-l 


i  r1  gn(?')  x  r1 


From  Eq.  3  the  second  term  in  Eq.  A-  1  becomes 


1  fl  gnKW  C0n  1 

— - = - 2*  Cvn  cos  v6  n  =  0,1,2 

Zt.  J-l  £  -  £'  2  2 


(A- 2) 


and  the  integral  expression  in  the  first  term  of  Eq.  A- 1  becomes 


J1  g„  (V)dV  = 


=  TrCr,„  +  - 


On  +  -  °li 
2 


(A- 3) 


The  last  term  in  Eq.  A- 1  is  then 


J  ^  gn(  =  ,)Un(r^d^'  =JQ  |C0uctg~+  v?i  '"vn  sin  v0' j  Un(q)sinG'  d0' 

=  JQ  jC0n+C0ncos9'  +  5l  Cvnsinv0'  sin  6'  J  Un{r))d6' 

=  |  ConUn(rj)d0'  +  j  Con cos  0' U^qJdO' 

Jo  Jo 

+  |  ^  Cvnsinv0'  sin0'  Un(q)d0'  (A-4) 

Jo  V  "  1 


From  the  identity 


2sinv0'  sin0  '  =  -cos(v  +  1)0'  +  cos  (v  -  1)0' 


(A- 5) 
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jT  gntt'JUnW  =  COn[£o(e)  +  fi(0)]+  £  C 


[*v-l<e>- We>l 


wn 


a0n 


an^)  =  — —  +  ^  apn  cos  p0 


(A-l  1 ) 


(A-12) 


From  Eq.  A- 7  and  A- 9 

2  r 7T  rn 


bV|jL  - 


.-ff 

ir  JO  J( 
2  /*tt  rT 

=  “  I 

tt  JO  JO 


cos  v0' Un(q)d9'  cos  |JL0  d0 


cos  vG'  cos  |jl0 Un(rj)d0 '  d0 
From  Eq,  A-2,  A-3,  A-ll,  and  A-12,  Eq.  A- 1  becomes 


(A- 13) 


a0n  n\  /  tt 


2  P 


+  0=1  apnCOS  P9  = —  -rrCon  +-  Cin 

4  2 


2  2 


—  i  C_  cos  v0  + 
v  =  1 


On 


2tt 


[fO(0)  +  fiO)] 


♦7&c«. 


fv-l(0)-fv+l(0) 


n  =  0,1,2  .  , .  (A- 14) 


a0n  y 

~  +  ptlapn«  = 


From  Eq.  A-9.  Eq.  A- 14  is  then 
n\  /  it 


\ 

C0n  1 

y  1 

)i 

2  2 

vtl  Cvn  cos  v0l 

^C0nj 

(b00 

2tt 

l  2 

b00  y  b10  y  \ 

—  +  fil  V  ':os  **  “  +  fil  bln  cos  & | 


^  £  Cvn 

2^V  =  1~T 

/bv+l,  0 


bv-1,0  x 


\  2 


+  p.  =  lbv-l,H  cosp-e 


+  I 

V-  = 


S  1  bv+l,  pi  cos 


(A- 15) 
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Equation  A- 15  can  be  written 

a0n  y  nW  77  \  /  ^0n  1  y 

—  +  p  =  l  apn  cos  P0  =  —  Kon  +  ~  Cln  +  I  —  ’  “  V=1  Cvn 


Cvn cos  V0 


XC0nlh00  |  b10^ 

\ 

“t* 

I  °vn 

bv-l,  C 

bv+l,  o\ 

2ir  \  2  2  i 

2tt 

V=1  2 

2 

2  1 

+  —  C0n  <b0jx  +  blp>  cosP6 
2  it  r 

^  v  0vn  v 

+  —  V  =  1  -  U=1  (bv-l,  M-  "bv+l,  |x>  cos  P0  (A"l6> 


Equation  A- 16  represents  the  general  boundary  condition  for  ring 
wings  in  series  form.  Let 


,?i  Cvn  cos  u0  =  p?i  Cpn  cosP° 


and  let  p  =  p.  in  Eq.  A-  16. 

Then 

n\-rr  n\ir 

a0n  +  2apncos  Pe  =  —  C0n  +  —  Cln  +  C0n  '  p?i  Cpncos  P0 

\  X  y 

+  —  C0n(b00  +  b10)  +  —  v=lCvn(bv-l,  0'bv+l,  0> 
2t r  4tt 


+  -  Con  (t>op  +  bip)  cos  pi 


+  —  v?i  Cvn  5i  (bv-l,  p  -by+1,  p)  cos  p 0  (A-17) 

2tt  p 

A  solution  to  the  general  wing  boundary  condition  can  now  be  rep¬ 
resented  by  equating  terms  in  Eq.  A-17  containing  cos  p9,  p  =  0,  1,  2  .  .  . 

For  the  case  p  =  0,  or  terms  independent  of  cos  p0 


20 


NAVWEPS  REPORT  8401 


irnX. 

X 

■tmX.  X 

a0n  ~  C0n 

1  + -  + 

—  (bio  +  b00^ 

+  Cln 

+  (b00  "  b20  ^ 

2 

2tr 

4  4ir 

+  ~  vtl  Cvn  (Vl,  0  *bv+l,  0)’  P  =  0  (A“18) 

4ir 


For  p  /  0  or  terms  dependent  on  cos  p8 


p?j  2apn  cos  p0  =  -  ^  Cpn  cos  P6  +  —  C0n  I  (bQp  +  blp)  cos  P0 


t r 


^  v  *^vn  y  , 

+  -  -  p=i  <bv-l(  p  ‘bv+l..  p)  cos  Pe  (A-!9) 

IT  2 

Isolating  the  term  fcr  v  =  p,  Eq.  A-19  becomes 

y  y  ^-C-On  y 

2apn  cos  p0  =  -  i“1  Cp.  cos  P0  + -  A  (bop  +  bjp)  cos  p8 


\  C 


pn 


+  _  (bp.1)p  -bp+1>p)  cosP0 


tr  2 


v  =  p 
X.  C 

+  -  V?1  —  p?i  (Vl,  P  "Vl,  y\  COS  p0  (A-2°) 

^V^p  2 


Expanding  Eq.  A- 20,  equating  corresponding  terms  and  dividing 
through  by  cos  p0,  p  =  1,  2,  ... 


Pn 


*^pn  ^On  ^ 

•—  +  ~  (b0p  +  blp)  +  ~  Cpn  (bp-l,  p  -bp+l,  p) 

2  2tt  4tt 


^  T 

+  v  =  i  0vn  (bv-l,  p  _bv+l,  p) 
477  v^p 


(A- 21 ) 


apn  ^pn 


1  \ 

*  T  +  ~  (bp-l.  p  'bp+F  p) 
2  4tt 


+  C0n —  (bj  +  bQ  ) 
Ztt  k  k 


\  Y 

+  —  v=l  Cvn  (bv-l,  p  "bv+l,  p)j  p  =  1-  2  • 
477  v,tp 


(A- 22) 


21 


NAVWEPS  REPORT  8401 


where  the  coefficients  bVp,  p  =  p,  are  defined  by  Eq.  A-7  and  A- 9. 

Equations  A- 18  and  A-22  represent  a  set  of  linear  algebraic  equa¬ 
tions  with  unknown  coefficients,  Cpn,  the  solutions  of  which  satisfy  the 
ring  wing  boundary  condition  in  Eq.  A-l.  A  solution  of  Cpn  is  available 
on  an  IBM  7090  computer  for  any  value  of  X.  and  any  component  of  local 
incidence  angle,  an(£)  represented  by  the  set  apn>  p  =  0,  1,  2,  ...  . 
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